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NOMENCLATURE 


beam 


L length of bean; 

B width of beani 

h height (thickness) of beam 

X,y,Z global co-ordinates 

u,v,w displacement along the X,Y,Z axes 

u displacement 

o 

^x'^y'^z strains 
stresses 

11 '^12' stiffness coefficients (Appendix - 1) 
16 


N ,M 

X X 


X 

e 


^ 1 * * '^4 


0l'02 

Ml-Mj 


0) 


F 

n 


force and moment aggregate 

density of composite 

local co-ordinate 

length of element 

shape fianctions 

natural frequency 

shear forces 

moments 

forcing frequency 
nodal force vector 

non-dimensional frequency parameter 



X 


9, Y, R 
? / n / C 
N^.N2,N2 

t . 

1 

V^,V2,V3 

i , j , k 

^L' \ 
^LT 

i 

L 

n 

S. . 

12 

t 

T 

u, v,w 

L 1 ^ 
[d] 

[b] 

|j1 

[s] 

[T] 

a 

e 

Y 


AXISYMMETRIC SHELLS 

global co-ordinates 
local co-ordinates of the ele’Dent 
shape functions 
thickness at node i 

unit vectors along !i,n directions 

unit vectors along 0, Y,R direction 

modulus of elasticity along L and T directions 

modulus of rigidity in L-T plane 

thickness of layer 

subscript used for nodes 1,2 and 3 

direction parallel to fibre 

total number of layers 

elements of stiffness matrix 

thickness of shell 

direction normal to fibre 

displacements at any point in the shell along 
Y, R directions 
transpose of the matrix 
stress-strain matrix 
strain-displacement matrix 
determinant of Jacobian matrix 
strain-stress matrix 
co-ordinate transfer matrix 

rotation of the normal to the shell in Y-R plane 
angle between iv and L 

displacement vector of the element 
normal strain 


e 



nornial stress 


shear stress 
shear strain 

Poisson's ratio measuring normal strain in T 
direction under uniaxial normal stress along L 
direction 

elemental stiffness matrix 
elemental mass matrix 
density of the composite 
natural frequency 
non-dimensional frequency 
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EPITOME 

A Computer Aided Design (CAD) software has been developed for 
the vibration analysis of Fibre Reinforced Plastic (FRP) laminated 
composite beams and a;< i symrriet r i c shells of revolution. 

For beams, laminate composite theory is used. The finite 
elerr>ent formulation for the vibration analysis of the laminated 
beaiTiS is presented. 

The formulation includes arbitrary number of bonded layers 
each of which may have different thicknesses and orientation of the 
elastic axes. 

In the case of laminated axisymmetric shells, isoparametr i c 
elements with four degrees of freedom per node including the normal 
rotation are used. The use of isoparametric element considers 
anisotropy across the thickness. The deformation of the normal to 
the mid-surface is included. Strain along the thickness is 
considered. The f orrriulation includes arbitrary number of bonded 
layers each of which may have different thicknesses, orientation of 
elastic axes. The formulation is compared for cylindrical shells 
with Flugge's theory and refined theory CET. Using this 
formulation, axisymmetric shells like conical shells, shells with 
Gaussian curvature and spherical caps can be analysed for 
axisymmetric modes of vibration. 

In the present software, analysis of free vibrations for 
natural frequencies and mode shapes and analysis of the response of 


sinusoidal forced vibrations are included. 


CHAPTER I 


INTRODUCTION 

Composite materials is the latest field in the development of 
the special purpose materials. These materials are replacing the 
conventional materials in every possible way. This is due to the 
fact that composite materi-ils can be fabricated according to the 
requirements of the design to meet the mechanical, electrical, 
thermal and other parameters. Present work deals with the field of 
fibre composites. The success of fibre composites results from the 
ability to make use of the outstanding strength, stiffness and low 
specific gravity of fibres such as glass, graphite, Kevlar etc. 
These outstanding mechanical properties are combined with the 
unique flexibility in the design capability and ease of 
f abr i cation . 

In the present work, a computer Aided Design (CAD) software 
is developed for vibration analysis of Fibre Reinforced Plastic 
(FRP) laminated composite beams and axisymmetric shells of 
revolution. Finite Element Method (FEM) is used. Analysis of free 
vibrations for natural frequencies and mode shapes and analysis of 
the response of sinusoidal forced vibrations are included in the 
software . 

LITERATURE SURVEY 

A detailed literature survey has been done on the following 


subj ects 
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Free Vibrations of Laminated Composite Axisymmetric Shells 

The composite laminated shells have found extensive 
applications in the aircraft and space craft industries. 
Considerable progress has been made during the last decade in 
applying the FEM for the analysis of shell structures. The FEIi 
applications in this area may be classified into vo types! 

1. Those in which triangular or quadrilateral elements have 
been used C103 . 


2. Those in which ring elements have been used E3,8,113. 

The first approach permits general treatment of the shell 
and necessitates the use of large number of degrees of freedom to 
simulate adequate inter-element compatibility. 

In the second approach certain known displacement 
distributions along the circumferential direction are chosen and as 
a result, the inter-element compatibility is easily achieved. But 
this introduces a limitation in the use of this element, for the 
basic shell problem, as it has to be reduced to a set of problems 
with suitable circumferential distributions. Generally a Fourier 
distribution of loads and displacerrients is considered. 


These 
frequencies 
frequencies , 
complicated. 


two types of elements are good for the analysis of 
in circumferential modes. For axisymmetric mode 
due to large number of elements, the problem becomes 
Till now, most of the work on the vibrations of 
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laminated composite shells is on circumferential mode frequency 
analysis. Application of FEM to axisymmetr i c mode frequency 
analysis is few. 

Analysis of circumferential mode frequencies is done in 
references C43 to C173. In ref. C43, C53 free vibration problem for 
circular and non-circular cross-ply laminated composite shells is 
studied. Shivakumar et.al C83 have used eight degrees of freedom 
per node for a typical conical laminated shell ring element. Wu 
et.al C9II used higher order curved shell elements. Bert et.al L\Q1 
studied the effect of various shell theories, thickness shear 
flexibility and bimodulus action on small amplitude vibration of 
thick circular shells laminate of bimodulus composite materials. 
Sheinman et. al C113 investigated that coupling between symmetric 
and antisymmetric inodes is a function of stacking combination S 
orientation. It’s influence decreases as the number of layers in the 
laminate increases. The influence of the coupling effect on natural 
frequencies decreases as the lowest frequency is approached. 
Sheinman et.al C12I1 applied a linear theory. Stavsky C133 used 
refined Love-type theory. Bert et.al C143 examined the non- 
syiTimetric modes of an arbitrarily layered general anisotropic 
finite length shell. Dong C153 used Donnel theory for free 
vibrations of laminated orthotropic cylindrical shells. Greenberg 
et.al C163 used Love-type theory for layered shells in which each 
layer can have fixed fibre orientation. Bhimaraddi C173 
concluded that the stress and displacement distribution across the 
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thickness of laminated composite shells is non-linear. This non- 
linearity increases with increase in the thickness and the material 
0 r t h o t r 0 py . 

The analysis of axisymmetric mode frequencies is done in 
reference C23 , C3I1 , CIST . 


In 

ref . 

L23 Sun 

et.al developed 

a refined 

theory for 

vibration 

of 

laminated 

composite thick 

shells. 

This theory 

includes 

the 

effects 

of the transverse 

shear 

deformation. 


transverse normal stress and strain, rotary inertia and higher 
order stiffness terms. Sankaranarayanan et.al C33 studied 
axisymmetric free vibrations of laminated conical shells with 
linear thickness variation in the meridional direction. Rayleigh- 
Ritz procedure is adopted for the analysis. Classical thin shell 
theory is used. Parametric studies are presented to illustrate 
the effects of geometric, material and coupling parameters and of 
the boundary conditions on the frequencies and mode shapes. 
Malhotra et.al II18I1 used a super parametri c curved element with 
four degrees of freedom per node including normal rotation for the 
vibration analysis of filament wound FRP shells of revolution. 
Effect of winding angle, metallic liner and material on 
frequencies is studied. 

Forced Vibrations of Laminated Composite Axisymmetric Shells 

Wu et.al C9I1 used higher order curved shell elements to 
investigate forced responses of thin laminated composite shells. 
The effect of geometric non-linearity was included in the 
determination of forced response using Newmark’s integration 
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scheme. Reddy CSSD presented an extension of the Sander's shell 
theory for doubly curved shells to shear deformation theory of 
laminated shells. The theory accounts for transverse shear strains 
and rotation about the normal to the shell mid-surface. Exact 
solutions for simply supported, doubly curved, cross-ply laminated 
shells under sinusoidal, uniformly distributed and concentrated 
point load at centre are presented. 

Free Vibrations of Laminated FRP Composite Beam 

Miller et.al C233 presented an analytical method for 
determination of resonant frequencies of generally orthotropic 
beams where flexural to torsional coupling is caused by the 

character of the elastic material constants. The analysis predicts 
that torsional as well as flexural resonant conditions may exist 

for a freely vibrating beam. Teoh et.al CE43 used a continuous 

model which includes shear and rotary interia. Abarcar 

et.al CE53 have formulated a discrete mathematical model based on a 
combination of Myklestad and Holzer's method of lumped mass to 
predict the natural frequencies and mode shapes of Fixed-free beams 
of general orthotropy. The model incorporates both transverse 
shear deformation and rotary interia effects. Teh et.al CEE3 
presented two finite element models for the prediction of natural 
frequencies of fixed-free beams of general orthotropy. The 

discrete models include the transverse shear deformation and the 

F 

; rotary interia effects. Teh et.al CE73 studied the effect of fibre 

I orientation on free vibrations of composite beams. The torsion- 

flexure coupling effect of generally orthotropic beams is dependent 
on reinforcing fibre orientation and mode order . At higher ranks 

i\ 

'Ht'' 

'SiF 



of vibration this coupling effect is principally contributed fay the 


twisting rrioment induced due to bending. The influence of fibre 
orientation on normal mode shape is more significant for small 
values of fibre orientation. 

The work cited above is on generally orthotropic beams. The 
survey indi cates that not much 1 i teratur e is available on the 
vibration of laminated composite beams. Singh Z263 has used finite 
difference method for laminated composite cantilevers for finding 
natural frequencies of beams with flaws and without flaws. He has 
also experimentally verified the results. 

PRESENT WORK 

The following problems have been attempted in the present 

work: 

1. Free vibrations of laminated composite beams. 

2. Sinusoidal forced vibrations of laminated composite beams. 

3. Free vibrations of lairiinated composite axisymmetric shells. 

4. Sinusoidal forced vibrations of laminated composite 

axisymmetric shells. 

Laminated FRP composite beams with following boundary 
conditions are analysed. 

(a) Simply supported 
(b> Fixed - Free 
(c) Fixed - Fixed 

Laminated composite theory is applied for formulating the 


finite element model. Galerkin method is used. 


The formulation 
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considers laminated beam with arbitrary number of layers and 
different thicknesses having arbitrary orientation of elastic aKes. 
If the laminate beam is constructed of stacking a number of 
orthotropic laminae in an arbitrary sequence of orientations the 
undesirable effects of coupling between bending and stretching will 


be present 

■ 

To avoid 

this , for 

each 

ply above the midplane 

an 

identical 

ply 

should 

be placed 

at 

equal distance below 

the 

midplane . 

This 

type of 

laminate is 

called symmetric laminate. 

Use 


of symmetric laminate will eliminate warping due to inplane loads. 

According to the Love— Ki r chhof f hypo thes i s , the normal to the 
middle plane remains normal to the deformed surface and moreover 
there is no stretching along these normals. This Love-Kirchhof f 
theory is applied to the present fo rmulat ion. It is also 
considered that the plane sections remain plane during the 
deformation. Thus bending-torsion coupling is not considered in 
the present formulation. 

In the formulation of laminated composite axisymmetric shells 
the effect of anisotropy along the thickness is considered. This 
anisotropy arises due to stacking of orthotropic laminae with 
different orientations in a laminate. Therefore Love-Ki r chhof f 
theory is not applicable specially in the case of thick shells. 

An attempt is made in present work to refine the mathematical 
model by discarding Love-Ki rchhoff theory and considering the 
following assumptions. 

(i) Normal to the middle surface does not remain straight after 


deformation. 
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(ii) Distribution of di s pilacements is non-linear and that of 
strain is linear, across the thickness. 

(iii) Normal to the middle surface does not remain normal after 
def ormat ion . 

(iv) Case Is Strains in the direction normal to the mid-surface 
are considered. 

Case IIs Strains in the direction normal to the mid-surface 
are not considered. 

In the finite element analysis shear deformation and 
coupling effects are included. An isoparametric element is used to 
include the effect of non-linearity of displacement across the 
thickness. Using this formulation a shell with arbitrary number of 
bonded layers having different thicknesses and orientation of 
elastic axes can be analysed. The element has four degrees of 
freedom per node including normal rotation. As the formulation is 
for a general second order parabolic curve, various axisymmetric 
shells like cylinder, cone, hyperboloid, cylinders with Gaussian 
curvature and spherical caps etc. can be analysed for axisymmetric 
mode vibration. 

Two types of boundary conditions are considered. 

(a) Simply supported without axial constraint. 

(b) Simply supported with axial constraint. 
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CHAPTER II 


ANALYSIS OF LAMINATED COMPOSITE BEAMS 
AND AXISYMI-IETRIC SHELLS 


LA-MINATED BEAM 


As shown in Fig. 2.1, the X, Y, Z axes are along 
length, width and thickness of the bea^. Let length, breadth 
and height of the beam be’L, B and h respectively. The 
deformation be u,v,w along X,Y,Z axes respectively. 

The bending deformation is in X-Z plane as shown in 
Fig. 2.2. Hence, v = 0 and u and w are constant across width. 
Stretching of the normal ABCD (Fig. 2.2) is negligible. Inplane 
strains are. 


e 


X 


9u 

8X 



Z 


2 

9 w 



(a) 


e 

y 


0, Y 



(b) (2.1) 


where , 


u = u 

o 


IE 

9x 


Z 


(c) 


displacement of the middle plane. 
In plane stresses are, 

'^x ^ ^11 ^x 

^y “ ^12 ^X 

'^xy “ ^16 ^x 


( 2 . 2 ) 






12 


where 0^2' are as given in APPENDIX-I, 


Aggregate of the forces and nonents for N-laysred laT-iinate 

N 


is (N^, M ) 




X 


K=1 Z 


a (1,Z) dZ 

X 


(2.3) 


k-1 


where Z^,, shown in Fig. 2.3. Substituting the 


value of CT we get, 


9u 

\ = ^-11 — - ^11 


a^w 

ax^ 


(a) 


(2.4) 


au 

M = B. 1 ^ 

X 11 3x 


- D 


11 


a^w 

3X^ 


(b) 


where, 


N 




-k 

^11 


[{Z^ - z^_^) , (z^ - z^_^) , 


.2 

■'k-1' 


<=> 


The equation for vibration of the composite 
beam is given by 


2 

3^M 


(x, t) 


P (x) 


3x 


a w(x.t> 
at2 


+ f (x, t) = 0 


(2.5) 


Substituting for M(x,t) we get 
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3x 


2 [P (x) 


where. 


3 w (x, t) T . . , . 
^ ^ _i + - (x) 

3x 


3 V (x, t) 


+ f (x, t ) = 0 

(a) 


( 2 . 6 ) 


B 


D (x) = 


11 

‘11 


- D 


11 


(b) 


Finite Element Analysis of the Beam 

A linear finite element with two nodes is taken 
with shape functions. 



e 

(2.7) 


where , 

x - local co-ordinate of the element, 
e - length of the element. 
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If w(x,t) is deflection at any point x and at time t 
as shown in Fip. 2.4 we have, 


w = LnJ (w ] , 


( 2 . 8 ) 


where , 


LNJ = L^l . Nj , N. , 


= deflection at node i 


Applying Galerkin method to equation (2.6a), we get, 


/ (Ln •j'^ J • D (x) - |J (x) L^J ^ 

r\ 1 


Ln'J'" [d(x) . ^3 ® - - LnJ 

8x 


^ ^ Cd(x) 


3x 


9x 


(2.9) 


/ [_nJ f(x,t) dx 
o 


Above equation can be written for an element as. 


( [k] - [M] ) {w^} = {F.} + 


+0- 

+M, 


( 2 . 10 ) 


The elemental stiffness and mass matrices are 
assembled and for free vibration case we get. 


^ W -■ M ) = 0 


( 2 . 11 ) 
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where, ["kI and [[mj are global stiffness and mas matrices. 
The equation (2.11) is solved using standard algorithm. 

The values of a. are the natural frequencies for the beam. 

If u is the forcing frequency of the sinusoidal 
load then we get, 

■; [k] - 5^ [m]) {q} = {F} (2.12) 


where, 

q - displacement at node 
F - nodal force vector 

Solution of equation (2.12) gives the maximum displacement 
at each node. 

The boundary conditions considered are as follows: 

(a) Simply supported 

w = 0 at X = 0 and x = L. (2.13) 

(b) Fixed-free 

w = 0 at x = 0 and w = 0 at x = 0 (.2.14) 

(c) Fixed-fixed 

w = 0 and w* = 0 at x = 0 and x = L respectively. 

(2.15) 


il 
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LAMINATED AXI SYMMETRIC SHELLS 

A general case of axisyrninetric shell is shown in 
in Fig. 2.5. Global co-ordinates 0, Y, R are shown with 
respect to the shell. The elenjents are shown on the cross- 
section of the shell with the respective nodes. This is a 
parabolic curved eleinent with three nodes per element and 
\ . riable thickness at each node. The local co-ordinates 
for each element are given by C ,n , s as shown in the 
figure. 

Considering a single element, the geometric shape 
functions are. 


~ (-11 + n ^) 

N^ = 1 - ri^ (2.16) 

N^ = (n + n^) 


The global and local co-ordinates for this element are 
related as follows: 


0 = 


{R} 


{Y } 

N { ■"}+! N. 


S V 


3i 


(a) 


(b) 


(2.17) 
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global co-ordinates at node i 

^ "h 

the thickness of the shell at the i node 

unit vector at node i in the direction 
normal to the middle surface 

1,2,3 (nodes) 

At any point in the shell along K , n, 5, V^, V 2 , are 
unit vectors. Let 1 , j, ic be unit vectors along 0, Y, R, 

Thus 

T* 

1 

(j . 1^ + ic . I^-) // I 

^ ^ I (3Y/3ti)^ + (3R/3Ti)2j 

{-j . 3R/3n + k.3Y/3n)/’^j (3Y/3 ti)^ + (sR/sn)^] 

(2.18) 

The unit vectors at node i are found by dropping the second 
term of right hand side of equation (2.17b). 

DISPLACEMENT FUNCTIONS 

At any node i the displacements are u^, v^, and 
the normal rotation in the Y-R plane as shown in Fig. 2.6. 



where. 


Yi,Ri 


t . 
1 


^3i 


1 ^ 4 - Element Number 

1,2,3 - Node Number 

9, Y, R - Global Co-ordinates 

r| - Local Co-ordinates 


Fig. 2.5 Global and local co-ordinate system 
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At any point on the shell the displacement u,v,w are 
given by, 


u 


V 


w 


i: N u. 
i 1 


IN. V. 

1 1 


IN. w. 
1 1 


+ IN. 

1 

+ I N. 

1 





a^. -f V 

j- 

a . + V 
1 


3i) 

3i^ 


where. 



1,2,3. 



components of unit vector V 2 at node 
along Y and R axes . 


Let, 


= L Vj_, w^, (a) 

For node, i = 1,2,3,. Therefore, 

« 2 ' “^3 

[u, v,wj'^ = C N 3 ■t (c) 

where [n] is matrix of shape functions. 


(2 


STRESS-STRAIN RELATIONS 


. 19) 


i 


. 20 ) 


The three dimensional strains with respect to global 
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axes are given by. 


^9 

1 

1 

1 

K + 1 

R R 

dU 

dO 

-J 

^Y 


3v/3Y 



^R 

= 

3w/9R 

1 


1 

1 

CD 


1 

■ i 3v 

R 96 

du 

9Y 


'^9R 


i ^ 

R 99 9R 

- 

u 

R 

1 


3w 9v 

9Y 9R 



L 

— 




( 2 . 21 ) 


The global derivatives are obtained as follows: 


9N, 

] 

99 


(a) 


9N. 
1 

9Y 


9R 

3? 


9N. 
1 

9 n 


/ 1 Jl 


Cb) 


( 2 . 22 ) 


9N. 

2 

9R 


9Y 


9N. 


3n 


- / IJ 


(c) 


where. 


J 


9Y 

9ri 


9R 

9n 


9Y 

3? 


IR 

3n 


(d) 
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The strains with respect to the local axes are 
obtained frori the strains along the global directions. 
Finally we get strains along the local axes in terTiis of 
displacements at nodal points as follows: 


Fe-e e Y- y 

- r r rH 




n?' 


B 1 { S 1- 


(2.23) 


Let 6 be the orientation of the fibres in a layer 
as shown in Fig, 2.7. Therefore, the transformation matrix 
is given as 


[t]- 


m^ n^ 0 -2mn 0 0 


m^ 0 2mn 0 0 


2 0 

mn -mn 0 m 0 0 


0 0 0 


0 0 


0 0 


m -n 


n m 


(2.24) 


where. 


m = sing , n = cosg 


The stress-strain relation with respect to elastic 
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axes is given by. 



(a) (2.25) 


where, 

L - Longitudinal axis of fibres. 
T - Transverse axis of fibres. 
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TDodulus of rigidity and Poisson's ratio. 

The stress-strain relation with respect to local 
axes are given by 



Sn 




]"=[t] [s]-" [t]^ e 


nc 




(2.26) 


The element stiffness matrix [k] is given by, 

[k] = fff [b]^ [t] [S]"^ [t]"^ [b] d (vol) 

11 

= f f [By [d] [b] 2TrR \J\ dn d? (2.27) 

-1 -1 

where. 


|j| = The Jacobian given by equation (2.22d) 

[d] = [T] [S]-l [T] 


As [b] matrix is different for each layer, the variable c in 
*tln 

the layer is changed to such that varies from 

-1 to 1 in the layer. Then, 


K 

5 = -1 * i [-he - V " ^ ^ ] 

\ 

d^ = “ ‘"'k 


(2.28) 
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where is the thickness of K^^'layer. With this substi- 
tution we get, 

r ^ 1 1 „ - h, 

W = s ( / / [b]^ [d]^ [b] 2TrR jjj ~ dn ) 

-1 -1 ’ 

(2.29) 


where n is total number of layers. 


Equation (2.29) is integrated numerically using a three 
point Gauss method. 

The elemental mass matrix is obtained as. 



Z f [n]T [p [n] d (vol) 

k=l vol 


n 

= Z 
k=l 



/ [n]'^ [pc^k Tn]2^R|J! 


^k 

— dn d^ 3 ^) 


(2.30) 


where, Fp J 

C"^ 


^ is the density of the composite in layer k. 


The element stiffness and mass matrices are assembled 
to get global stiffness and mass matrices. They are used to 
solve the generalised eigenvalue problem. 


[k] {6} = 0)^ [m] {5} 


(2.31) 


The values of oi give natural frequencies and the values 
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of {5}give mode shaoes. Equation (2.31) is solved 
using standard eigen value solution subroutine. 

The sinusoidal forced vibration resoonse is 
obtained as given by equation (2.12). 

The two types of boundary conditions considered 
for the axisymmetric shells are as follows: 

(a) Simply supported without axial constraint 
u=0, w=0, a =0 and v?^0 at both ends. 

These boundary conditions can be approximated by 
attaching thin, flat, circular diaphragms to the ends 
of the cylinder. The diaphragms have considerable 
stiffness in their own plane, thereby restraining 
the u and w components of the shell displacement at 
their mutual boundaries. However, the diaphragms, 
as a result of their thinness, generate negligible 
bending moment and longitudinal force in the shell as 
the shell deforms and do not restrain axial deformation, 

V. 


Simply supported with axial constraint 
u=0, v=0, w=0 and a=0 at both ends. 


COMPARSION WITH OTHER HETHODS OF ANALYSES: 


Laminated Beams 

Most of the work is done on the orthotropic beams but on 
anisotropic laminated FRP beams, work is rather scarce. Singh 

C26I} used laminate composite theory with finite difference method. 

In the presen't, work, the laminate theory is used with the 
following assumptions. 

1. The normal to the middle surface remains normal to the 
surface after deformation. 

2. The stretching along this normal is negligible and hence 
taken to be zero. This means that there is no strain along 
the thickness direction. 

3. Bending deformation is in a plane and hence displacement 
along width is zero and displacement along thickness and 
length are constant across the width. 

4. Plane sections remain plane during the deformation. 

The present model does not consider shear deformation 
and bending-torsion coupling. 

Axisymmetric Shells 

Most of the work on vibrations of laminated fibre reinforced 
composite shells of revolution considers the circumferential mode 
frequency analysis. Therefore either triangular or quadrilateral 
elements have been used. This led to consideration of large number 
of degrees of freedom per node to simulate adequate inter-element 
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compiat i bi 1 i ty . If these elements are used for axisymmetric mode 
frequency analysis, the size of stiffness matrix becomes very 
large . 

Conical and cylindrical elements were used with Fourier 
distribution of loads and dis placements along the circumferential 
direction. These elements are useful for getting both types of 
frequencies for cylindrical and conical shells. But these elements 
involve approximations when used for general axisymmetric shells. 

The use of parabolic curved shell element in the present 
analysis is suitable for general axisymmetric shells. Some of its 
features ares 

1. The approximation for the surface is not required. 

2. The axisymmetric thickness variation along the length of 
shell is permitted. 

3. Variation of density and elastic coefficients for each 
lamina in the laminate is permitted . Due to the use of an 
isoparametric element following additional features arise. 

4. Distribution of displacement is non-linear across the 
thickness and that of strain is linear. 

5. Due to lamination of orthotropic layers, anisotropic 
properties are present across the thickness. Therefore, a 
refinement is made from other theories in assuming that the 
normal to the middle surface after deformation does not 
remain straight. 

6. Strain in the direction normal to the mid-surface is 


considered. 
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CHAPTER III 

RESULTS AND DISCUSSION 

The software developed has five different modules. The 
execution of the main program "PROJECT" will lead to the MENU. In 
this MENU, all modules are connected as shown in Fig. 3.1. 
Selection of a particular module can be done with the help of a 
cursor. Functionally the modules can be catagorised in the 
following way. 

(I) Laminated Beams 

There are two modules for the analysis of laminated 

beams . 

(a) Vibration analysis of laminated beams 

The vibration analysis can be done for two cases; 
(i) Analysis of Natural frequencies and mode shapes. 

<ii) Analysis of response for the sinusoidal forced 
vibration . 

Module "FVBEAM" is for the vibration analysis of the 
laminated composite beams. It works in the following different 
phases . 

(1) The input data is collected from the user on the terminal. 

This data is also written in the file "DAT-INP". 

(E) The processing of data is done and the coefficients A-^-|, 
from equation (2.4) and D(x) from equation (E.6 b) 
are calculated. 

<3) Elemental mass and stiffness matrices are calculated. These 
are then assembled according to the boundry conditions to 
form global mass and stiffness matrices. The equation 
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(2.11) is solved using NAG subroutine (F02AEF). The results 
i.e. natural frequencies and modal vectors are stored in 
the file "DAT-OUT". 

(4) For getting the sinusoidal forced vibration response, first 
a file "DATA: INF" is to be created in which the forcing 
frequency is written first and then the nodal values of 
force are written. The response vectors are written in file 
"DAT-OUT". 

(5) Plots of mode shapes and response are then drawn on the 
workstation. 

(b) Static load analysis of laminated beams 

"STBEAM" is for static load analysis. The execution of 
this program will give the strains and stresses in each layer at a 
particular section of the beam. Different types of loading 
conditions are provided. This program is based on the exact 
analysis. Execution of "STBEAM" creates "SAT;INP" and "SATsOUT" 
files. 

(II) Ax i symmetric Shells 

There are three modules for the analysis of laminated 
axisymmetric shells. 

(a) Execution of "INCUR" will lead to an interactive 
session with the user . User has to feed data on the 
terminal. This data consists of following items: 

(i) Element node point co-ordinates in Y-R plane as 
shown in Fig. 2.5. Thickness at each node. 

(ii) Material properties. 
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(iii) Number of laminae, orientation and thickness of 
the laminae. 

(iv) Forcing frequency cf the sifiuscidal force and the 
nodal force values. These values will depend on the 
type of loading e.g. concentrated load, unif ormely 
distributed load etc. 

(b) Module "CURVE" is for the vibration analysis of 
laminated composite axisymmetric shells. 

(c) Module "ZETACUR" is for the vibration analysis of 
laminated comiposite axi symimet r i c shells in which the 
strain along the direction normal to the middle plane 
is also considered. 

After execution of any of these two modules i.e. "CURVE" 

and "ZETACUR" following operations are performed; 

(1) Elemental stiffness and mass matrices are calculated 
using equations (2.29) and (2.30). 

(2) Natural frequencies and modal vectors are obtained 
using the equation (2.31). Natural frequencies are written 
in file "CURVE50UT". Modal vectors are written in file 
"COUT". 

(3) The response is written in file "CURVE:0UT". 

(4) The plots of mode shapes are drawn on the workstation. 

(5) The plots of response in the direction 0,Y,R are drawn 

separately on the workstation. 

For the detail knowledge about the program please refer to 




Appendix II. 
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RESULTS FOR COMPOSITE BEAMS AND SHELLS 
Laminated Composite Beam 

The program "FVBEAM" was checked for the non-dimensional 
frequency parameter slj , for the isotropic case. The non- 
dimensional frequency parameter Qj is related to the natural 
frequency as ,• 

/ L^ (El/f) j = 1, 2, 3 (modes) 

In tables 3.1, 3.2, 3.3 a comparison is given for the results 
obtained from "FVBEAM" with exact results from reference C213. The 
convergence of the eigen values with increase in the number of 
finite elerrients to the exact value can be seen from the tables. 

TABLE 3.1 Comparison of eigenvalues ( S'ij ) for Simply Supported Beam 


Mode 2-Elements 4-Elements 8-Elerrients 10-Elements Exact 


Numbe r 
j 

1 

9.9086 

9.8730 

9.8822 

9.8208 

9.8696 

2 

43.8178 

39.6347 

39 . 4930 

39.4830 

39.4784 

3 

110.1396 

90.4495 

88.9408 

88.8715 

88.8264 

4 

200.7984 

175.2713 

158.5369 

158.1729 

157.9137 

5 


— 

249.0268 

247.7139 

246.7401 
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TABLE 3.2 Comparison of the eigenvalues ( i'ij) for Fixed-Free beam 


Hode 

Number 

j 

2“E1 ement 5 

4-El ements 

8-Elements 

10-El ements 

Exact 

1 

3.5176 

3.5163 

3.5670 

3 « 0995- 

3.5008 

2 

22.2214 

22 . 0605 

22.0388 

22 . "304 

22 . 0345 

3 

75. 1571 

62.1751 

61.7354 

61.7198 

6 1 » 697Z 

4 

216. 1360 

122.6577 

121 . 1729 

121.0151 

120.8990 

5 

— 

— 

201.0170 

200.3607 

199.8595 

TABLE 

3.3 Comparison of the eigen 

values ( Oj ) for 

Fixe d— F i x e d 

beam . 

Mode 

4-Elements 

6-Elements 

10~“E1 ements 

Exact 


Number 

j 






1 

22.4029 

22.3768 

22.3789 

22.3733 


2 

62.2433 

61 .7122 

61.6856 

61.6728 


3 

123.4654 

121.1911 

121 .0237 

120.9034 


4 

233 . 6242 

201 . 1068 

200.3897 

199.8594 


5 

366.3797 

302.4677 

300.2653 

298.5555 



For high strength Graphite Epoxy Composite material following 
input data is taken. 


El 

= 1.2 

X 

10^ 

N/mm^ 

'^^LT = 

0.27 

Et 

= 1.0 

K 

10^ 

N/mm^ 

L 

200 trim 

GlT 

= 5.5 

X 

10^^ 

N/mm^ 

B 

30 mm 

f 

= 1.7 

X 

10-*^ 

2 -4 

N Se c mm 

h 

5 mm 
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For the laminated beams, C0,90D5 and C45,-45D5 composites are 
considered. Table 3.4, 3.5, 3.6 shows the values for the natural 
Frequencies in rad. /sec. For three types of boundary conditions. 

(a) Simply supported beam 

(b) Fixed-Free beam 

(c) Fixed-Fixed beam 

The natural frequencies are plotted Vg the rriode number for 
CO, 903s C45,-453s composites in Fig. 3.2, 3.3, 3.4. It is 

observed that the frequencies for the CO, 903s composite beam are 
higher for each miode than that of C45,-453s composite bearri. 
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TABLE 3.4 Natural Frequencies For Laminated Composite Beams. 
(No. of Finite Elements = 10) 


SIMPLY SUPPORTED BEAM 


Mode 

1145, -45Ds 

C0,90Ds 

1 

469.26 

518.82 

H 

1882.96 

2064.13 

3 

4237.79 

4643.07 

4 

7542.49 

8263.49 

5 

11811.77 

12941.03 

6 

17076.76 

18709.32 

7 

23387.19 

25263.18 

8 

30813.37 

33759.29 

9 

39409 . 50 

43177.22 

10 

52234 . 66 

57228.49 

11 

61856.63 

67770 . 48 

12 

75400.48 

82609.19 
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TABLE 3.5 Natural Frequencies for Laminated Composite Beams. 


(No. of Finite Elements = 10) 


MODE 

FIXED- FREE BEAM 

C45,-453s 

CO, 901s 

1 

167.61 

171.86 

2 

1050.00 

1150.20 

3 

2942.85 

3223.99 

4 

5770.67 

6321.99 

5 

9554.02 

10467.31 

6 

14313.03 

15681.27 

7 

20081.92 

22001.72 

8 

26904. 18 

29476 . 30 

9 

34786.72 

38112.52 

10 

43244.93 

47379.28 

11 

57557.07 

63059 . 78 

12 

69536. 14 

76184. 16 
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TABLE 3 


'li' 



6 Natural Frequencies for Larriinated Composite Beams, 


(No. cf Finite Elements = 105 


MODE 

FIXED-FIXED 

C4D,-45>5 

BEAM 

CO, 900s 

1 

1066.88 

1 1 68 . 78 

z 

PQ4 'i 

3222.54 

3 

5770.75 

6322.24 

4 

9555.19 

10468.39 

5 

14371.62 

15686.46 

6 

20096.54 

22017.84 

7 

26943.22 

29519.21 

8 

34875.34 

38209.62 

9 

43384.16 

47531.88 

10 

57886 . 75 

63421.02 

11 

70232.94 

76947.56 

12 

85485 . 66 

93658.50 
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Axisytrifiietric Shells of Revolution 

In the program "CURVE", the equation ( 2 . 3 :) is solved using 
NAG subroutine (F02AEF) for getting the natural frequencies and 


modal vectors. This subroutine arranges frequenices 


t f'^ e 


ascending order. It is geriera^.l obsarvetic' (alsc ref .Ei ) that the 


natural frequenices of shells do not fall in an ascend 


H i 


Fi Q 


series 


e a c fi 


VMith increasing values of the madal index. Therefore 
frequency has to be checked for the modal shape. This will require 
lot of executions of "CURVE" module. The importance of the 
graphical plotting of the mode shapes is strongly felt while noting 
the natural frequencies. The numerical results are obtained from 
the program "CURVE". For high modules Grapihite Epoxy composite, 
following data is 


given as input, 


-L 


= 20 X 10^ psi 

Gtt’= 0.5 X 10*^ 

Et 


= 1 X 10^ psi 

= 0.25 

GlT 


= 0.6 X 10^ psi 

^j-r = 0.25 


The 

Tion*“dimens i onal frequencies are given by y 


Q = 

^h/An'T'^/EL where 

A-y\ = n w h/L 


n 

= Wave number (half) 



L 

= Length of shell 



h 

= thickness of cylinder 



f 

= density of composite 



U3 

= natural frequency (rad 

/sec) 


R 

= radius of cylinder 



ps: 


I 
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For cylindrical 

shell. 

the 

plots 

f c r riori-difTieris 

i o n a. 1 

f r e 

quency Q 

are 

shown 

in Fig. 

wpOf 'Uj*/, 

3 . 8 y . 

For 

two cases of R/h = 

10 and 

R/h = 

3 the 

results are taken 

f G r 


CO, 90^5 cross-ply laminated shell and C45, -4501c angle-ply laminated 
shell of the composite material. 

Comparing the plots with those of reference C23, following 
observations can be made. 

'i- Fig. 3.5 to Fig. 3.8, the plots are between Q the non- 

dimensional frequency and /tr . 

As ‘'^n/TF > 0, h — > 0 but in the expression fof 0, 

52 ='^L/tt nJ" ^/Ei_. This means that as ^n/xr — > 0, 1'2 can have any 

value which is finite. 

2. In the present case the variation of the natural frequency 
with the mode number n is approximately linear as observed from 
Fig. 3.9 and 3.10. 

3. For a particular value of -^n/Tr = constant, as the mode 

number n increases the ratio (h/L) decreases in the same 
proportion i.e.1/n ;< (h/L>. This can be accomplished in two 
ways. One way is by decreasing the shell thickness h and keeping 
L constant. The other way is to keep the shell length L 

increasing and keeping h constant. In both ways the shell 
approaches to thin and long shell. 

Presence of four degrees of freedom at each node leads to 
large stiffness coefficients for relative displacements along an 
edge corresponding to the shell thickness. This presents numerical 
problems and leads to the ill - conditioned equations when shell 
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til i cknesses become small compared with the other dimensi ons in the 
element. Therefore in the present case, for thin, long shells the 
results are not relisible. 

In the case of thick shells with R/h = 3 the results are 
cotTiparable with those of reference ZS3. For ■^n/tr >0.5. In case of 
R/h = 10, the results are corri parable with reference LZl for 

Xn/TT > 0.2. 

4. The results obtained from Flugge's theory are higher 

than that of refined theory of ref. 023. Flugge’s theory considers 
only inplane strains. The refined theory 023 considers the 
transverse normal strain and the effect of transverse shear 
deformation is included by introducing an approximate shear 
correction factor. The present mathematical model considers all 
three dimensional strains. The effect of anisotropy across the 
thickness, the non-linear distribution of displacement across 
the thickness and the linear distribution of strain across the 
thickness is considered. All these considerations lead to a 
refined model. Therefore, the non-dimensional frequencies 
obtained in the present case are smaller in magnitude. 

From the analysis of the results it was observed that 
the natural frequencies change with the orientation of the fibres. 
A close observation at the pilots of Fig. 3.11 leads to the 
following point. 

The frequencies for the C0°3 composite are higher than that 
for the other orientation angles up to C45°3. At the 
orientation of [:45°3 the frequencies are minimutn- As the 
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Fig. 3.5. Nondimensional Natural Frequency as a 
Function of h/l. 













1400 



CO 1200 



Fig. 3.9 


Natural Frequency vs. Mode Number . 
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orientation angle is further increased from to C90^3 the 

frequencies are beconarig larger, specially for liigher modes. 

The reason for this kind of behaviou*' is obvicus when we 
consider the fact that the stiffness of the Fibres along 
longitudinal direction is higher than that in transverse 
direction. Therefore the composite shell with fibre orientation 
of C90°I3 is stiffer along the circumferential direction. This 

does not allow the shell to expand and contract according to the 
ax i symirriet r i c rriode shapes. And for the higher modes this 
resistance is substantial. 

To confirm this, for 11011° and [1901]° fibre orientation the 
transverse modulus of the fibres is increased from 

£[_/£■}- = 20 to E;_ /E j = 10. The results are surrimar i sed in TABLE 
3.7. From the study of this table, it can be observed that as 
transverse modulus Ej increases the frequency also increases. The 
frequencies for C90°3 fibre orientation are higher than that 
for C0°D fibre orientation. 

While executing the prograrii "CURVE" for ax i syrrimet r i c 
shells following observations can be made. 

1. The lowest frequency is not the first mode frequency. 
This observation confirms with that of ref.C213 < page 292). 

2 . It is observed sometimes that for a particular mode 
there are two frequencies. One frequency has very low 
amplitude of the modal vector while for the other frequency 
the ajTiplitude is higher. This is a characteristic of 
composite shells. Due to the bending-stretching and bending- 
torsion coupling these frequencies arise. 
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3. It is observed that bending-s t r et ch i ng touplirg --f-fetts 
the mode shapes in the esse of the bsinds'-y L:;:diticr. i 



s i m p 1 y 

supported 

shell 

w i t a i a 1 c d n r. 

* y. ^ y 4- 

i O. i. ’ I V « 

TABLF 

O “7 

Effect of 

chang i 

r in Ej . ^4atu^al 

F r e q uen c i es v^) rad/sec 

MODE 


FIBRE 

0RIENTATI0^4 !:0°3 



El 

= 20 


"l = -is 

El = '3 


T 




E,j 

•1 

49 

.62 


265 . 1 -1 

723.07 

2 

92 

.80 


305.2-1 

740 . 92 

3 

•136 

.05 


334.82 

764.01 


MODE 

FIBRE 

£[_ = 20 
Et 

ORIENTATION C90°3 

El = '‘S 

E? 

Er_ = 10 

Ef 

-1 

■1828.65 

2280 . 00 

3420.13 

2 

1831.65 

2282 . 00 

3422.86 

3 

1836.67 

2287.01 

3427.39 





55 


The present FEI^I formulation was reduced to the 
assumptions of Flugge’s theory by considering in-plane stresse 
and strains only. These results are Tabulated below. 

'lable 3.8 Values of Q (No. of elements =6) 


Kn 

n 

Flugges 

Theory 

^1 

Present 

FHyi 

s 


0.2 

0.1050 

0.1089 


0.4 

0.2036 

0.2073 

[0,90 Ig 

0.6 

0.3080 

0.3099 

h 1 

0.8 

0.4007 

0.4089 

R = To 

1.0 

0.5117 

0.5221 


1 .2 

0.6220 

0.6343 


0.2 

0.2555 

0.2918 


0.4 

0.5075 

0.5191 

[45,-45 ]g 

0.6 

0.7511 

0.7633 

It 

o 

0.8 

1 .0193 

1 .0333 

1 .0 

1 .2508 

1 .3233 


1.2 

1 .5009 

1 .5636 
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Flugges 

Present 


71 

Theory 

FEK 


0.2 

0.2564 

0.4117 


0,4 

0.5111 

0.6197 

_ 3 

0 .6 

0.7667 

0.8106 

[0,90 Ig 

0.8 

1.0135 

1 .2313 


1 .0 

1 .3333 

1 .4186 



0.2 

0.1511 

0.2189 


0.4 

0.2525 

0.3331 

& = 3 

K 

0 .6 

0.5010 

0.6333 

[45,-46 ]g 

0.8 

1 .0900 

0.8778 


1 .0 

1 .4500 

1 .1998 
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CHAPTER IV 
CONCLUSION 

Vibration analysis is done for the laniinated composite 
ax i symmetric shells and beams. 

Axisymmetric Shells 

An isoparametric curved shell element with non-linear 
distribution of displacement and linear d : st r i fcu* i or. of strain 
along the thickness is used. All three dimensional strains are 
cons i do red . 

From the analysis of the results following conclusions can 
be drawn. 

1. Variation of frequencies with mode number r is found to be 
1 inear . 

2. The frequencies for C9C°D orientation are found to be higher 

than that for orientation. This behaviour is more 

pronounced in the case of higher modes. 

3. The importance of the above-mentioned factors has been 
demonstrated by comparing solutions froiri the shell theory 
for contemporary engineering laminates to solultlions 
obtained from Flugge’s theory and refined theory of 
ref .C23. 

4. The effect of consideration of strain along thickness on the 
results in not substantial. 


Beams 


Laminate composite theory is used for 


the 


vibration 
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analysis of laminated composite beam. Three types of faoundry' 
conditions are considered namely; 

a. Simply supported beam 

b. Fixed-free beam 

c. Fixed- fixed beam 

The frequencies for CO, 9013s laminated beam are found to be 
higher than that of CO, 453s laminated beam in all three types 
bound ry conditions. 

Scope for Future Work: 

In the present work, in case of beams laminated composite 
theory is atpplied. In the case of laminated axi syrrimiet r i c shells, 
stress is given on the use of iso paramet r i c element. The following 
scope for future work exists: 

1. Laminated Beams 

a) Inclusion of anisotropy along the thickness. 

b) Inclusion of rotary shells 

E. Laminated ax i sy'mmet r i c Shells 

a) Analysis of results in the case of sinusoidal forced 
vibrations . 

b) Analysis of results in the case of conical shells, 
shells with Gaussian curvatures and other 
axi symmetric shells, spherical caps etc. 
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APPENDIX 
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An orthotropic lamina with principal axes orie 
at angle © with the reference axes is shown in figure. 
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where [t] is transformation matrix given by 
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11 - ^11 cos^e + O 22 sin'^G + 2{0^2 + 20g5) sin^G cos^G 

22 “ ^11 sin^G + ©22 cos^G + 2 ^0^2'^20gg) sin^G cos^G 

12 ~ ^^ll'*’®22~^^66^ sin^G cos^G + 0^2 (cos'^G + sin"*e) 

56 ~ ^^ll'‘'®22 ■ ^^12 “ ^^66^ sin^G cos^© + Ogg (sin'^O+cos^e) 

^16 ^ (021~^12"^®66^ ~ ^^22“^12“^^66^ sin^G 

^26 ^^11"®12"^®66^ sin^G - {022-0^2~^^66 ^ 
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The present software PROJECT: FORT is based 
on the ND-560 system. The FORTRAK-77 Language is used for 
the programming. The software uses available M4G and GKS 
library subroutines. 

The various modules connected to the main 
program PROJECT: FORT are as shown in Fig, 3,1. These are : 

(1) INCUR : 

This subroutine while executing, collects the 
data from the user which is used as input data for the 
another module CURVE, The following data is required : 

1 . The number of elements for the axisymmetric shell 

2, The geometry of the element described by the nodal 
points in the Y-R plane as shown in Fig. 2,5, 

3, Thickness of the element at each node 

4. Material properties as follows : 

Longitudinal modulus of the composite 
Transverse modulus of the composite 
Shear modulus in L-T plane 
Shear modulus in T-T plane 
Poisson's ratio in L-T direction 
Poisson's ratio in T-T direction 


% 

\t 

%T 

’'tt 
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5 . 

6 . 

7. 

8 . 

9 . 


10 . 


:niser cf layers in the ccr.aos 


aniCKnesi 


i-.ngle of orientation for each layer and the 
thickness . 

Density of composite. 

Boundary conditions for shell . 


Forcing frequency in case of forced sinusoidal 
vibration . 

The nodal force vector. 


The file CURVE : INP is written by the execution 
of INCUR . 


(II) CURVE : 

This subroutine is for the vibration analysis of 
laminated Fibre Reinforced Plastic composite axi -symmetric 
shells , The working of subroutine is as follows : 

(1) [s] matrix is calculated according to equation (2.25) 

(2) [s]”"' matrix is calculated using (F01AAF) NA.G 

subroutine . 

(3) Then according to equation (2.29) and (2.30) the 
elemental stiffness and mass matrices are calculated. 
The subroutine SHAPE or SHAP is called according 

to the coiKiition of 

1 . The strain along thickness is neglected - 
Module CURVE 

2. The strain along thickness is considered - 
Module 7.ETACUR . 
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The suhroutirie; 


n.-.rr. o: 


celcul' 


matrix 


L J (."j 


using the shane functions 


(^) The subroutine BOUND assigns the boundary 

conditions for the Global mass and stiffness 
matrices , It also decides their size by following 
expression. 


NX= (2x HOE +1) X 4 - KBC 

where 

NX = total number of degrees of freedom 

NOE = number of elements 

KBC = known boundary conditions 

(5) The subroutine ASSEMB assembles the elemental 

stiffness and mass matrices into the Global 
stiffness and mass matrices according to the boundary 
conditions . 


(6) The subroutine (F02AEF) is called from NAG library. 
This subroutine solves equation (2,31), The output 
is stored in files CURVE : OUT , Eigen values 
and vectors are obtained as output, 

(7) 


After selection of particular frequency by user 
the vectors are saperated for this frequency. This 
is done by using boundary conditions. 
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(£) 


(9) 


These vectors give "the disrlscerero at esoi. r.oce , 
These give mode shapes which are then drawn using 
GKS subroutines . 

Subroutine FORCE is for the sinusoidal forced 
vibration response. This subroutine solves the 
equation (2,12) using the subroutine (FOlAAF) from 
NA.G library. 


The response obtained is then processed 
and drawn on the screen using GKS, 

The module CURVE was executed for the input 
data given below : 



20x 10^ 

psi 

It 

lx 10^ 

psi 

^LT" 

0.6x 10 

^psi 


0.5x10 ^ 

psi 


^LT “ ’^TT 

R = 0.5 , 1= 10 

to. 903^ 


The mode shapes obtained are shown in Fig. 1 . The 

-3 

normalization factor is 2,73657 x 10 , 

For the sinusoidal forced vibration problem the 
following input is given: 

Concentrated force at node 1 of each element except the 

11 

first element == 2x 10 lb. 
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Frequency of forced vibretion = 2x 10^ rad /sec. 

The force is in positive R-direction. 

The response in 9 and Y direction is negligible (10”^ "^--10**^”^ 
The response in the R-direction is as shown in Fig, 2, 

(III) ZETACUR : 

This module works like CURVE except that it 
considers the strain along the thickness direction 
of the shell. 

(IV) FV BEAM : 

The execution of module FV BEAM leads to calling 
the subroutine FV BEAM. This subroutine is for the vibration 
analysis of laminated composite beams. The working of this 
subroutine is as follows : 

(1) The input data is collected from the user. This 
data consists of the following items : 

(i) The length, breadth, height of beam. 

(ii) The material properties of the composite 
laminate . 

(iii) Number of finite elements. 

(iv) Boundary conditions . 

(2) Subroutine BOUNDER is called for assigning the 
boundary conditions. The size of the global matrix 
is decided by the following expression : 


O- Cl ^7 
0 ■0‘^T 


o-oa.1 



1 I \ j I 

® © ® © © g 


Fiq 1 


R^spon^e 1« K- -C'S-ecitiom 
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I’K = (FELENO +1 )x2 - known bounoary conditions . 
where, NX = total degrees of freedom 

FELENO = number of finite elements. 

(3) Subroutine ASSMBL is^ called for assembling the 

elemental mass and stiffness matrices into global 
mass and stiffness matrices, 

(4) Subroutine F02AEF is called to solve the eigen 

value problem from NAG library. 

(5) Subroutine FORSIN is called for the sinusoidal 

forced vibration problem which uses FOIAAF 
subroutine from NAG library to solve equation (2.12). 

(V) ST BEAM : 

This module calls subroutine ST BEAM . This 
subroutine takes the input data as : 

(i) Geometrical dimensions of beam 

(ii) Material properties of the composite laminate 

(iii) Loading types 

(iv) Boundary conditions 

(v) Point at which the stresses and strains are 
to be found out, 

(vi) Points at which the bending moment and shear 
forces are to be calculated. 

The various subroutines are for the different 
loading conditions and boundary conditions. Exact analysis 
is done. 
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iu4238 

^ ^ lf?ok Is to be returned on the 

date last stamped. 
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